Spectrum of magnetohydrodynamic turbulence 
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We propose a phenomenological theory of strong incompressible magnetohydrodynamic turbulence 
in the presence of a strong large-scale external magnetic field. We argue that in the inertial range of 
scales, magnetic-field and velocity-field fluctuations tend to align the directions of their polarizations. 
However, the perfect alignment cannot be reached, it is precluded by the presence of a constant 
energy flux over scales. As a consequence, the directions of fluid and magnetic-field fluctuations at 
each scale A become effectively aligned within the angle <j>\ oc A 1//4 , which leads to scale-dependent 
depletion of nonlinear interaction and to the field-perpendicular energy spectrum E(k±) oc k ± 3 ^ 2 . 
Our results may be universal, i.e., independent of the external magnetic field, since small-scale 
fluctuations locally experience a strong field produced by large-scale eddies. 

PACS numbers: 95.30.Qd, 52.30.Cv 



1. Introduction. — Magnetohydrodynamic (MHD) tur- 
bulence is pervasive in astrophysical systems, where 
ranges of scales available for plasma fluctuations span 
many orders of magnitude, and the fluctuations com- 
monly possess power-law spectral distributions, e.g., 0. 
For example, the spectrum and structure of MHD fluc- 
tuations are relevant for the physics of solar wind, inter- 
stellar scintillation, cosmic-ray propagation in galaxies, 
and heat conduction in cooling flows in galaxy clusters. 

The spectrum of MHD turbulence was first addressed 
by Iroshnikov and Kraichnan who proposed the 
physical framework for describing the turbulent energy 
cascade mediated by a guiding magnetic field. How- 
ever, recent numerical and analytic works have chal- 
lenged these standard results and revived substantial 
interest to the fundamentals of strong MHD turbu- 
lence HUES, SHEIllElEl To formulate 
the problem and to set the notation, we first describe the 
Irosnikov-Kraichnan and Goldreich-Sridhar Q the- 
ories, and point out some discrepancies of these theories 
with recent high-resolution numerical findings. Then, we 
propose a new model for MHD turbulence, which is free 
of such discrepancies, and which explains the results of 
numerical simulations [5 I^L HoL Hi) . 

Consider a conducting fluid stirred by a random force 
with the correlation length Ao. The system size is larger 
than Ao, and viscosity and resistivity of the fluid are very 
small. The goal is to find the stationary energy spectrum 
of the resulting turbulent fluctuations in the inertial in- 
terval of scales, A <C A . Let us split the magnetic field 
into two parts, B(x,t) = B 4- b(x,t), where B = (B) 
is the system-size averaged magnetic field, and h(x, t) is 
the fluctuating part. The MHD equations describing the 
evolution of the magnetic field and of the fluid-velocity 
field v(x,i) can be represented in the so-called Elsasser 
variables, z = v b, and w = v + b: 



Stw - (V A • V)w + (z • V) 



-VP, 



(2) 



where = Bo/V^^P is the Alfven velocity, p is the 
fluid density, P is the pressure that is determined from 
the incompressibility condition, V-z = 0orV-w = 0, 
and we omit the terms representing large-scale forcing 
and small viscosity and resistivity. 

To present the standard arguments of Iroshnikov and 
Kraichnan, let us note that owing to the symmetric form 
of system two classes of exact solutions exist. For 

w = 0, any function z = g(r — V^t) is the solution of 
the system; analogously, for z = 0, the solution is given 
by an arbitrary function w = h(r + V^f). From the 
form of the nonlinear terms in system Q |2J , one observes 
that Alfven-wave packets, or "eddies," propagating in the 
same direction along Bo do not interact. One has there- 
fore to investigate interactions or eddies propagating in 
opposite directions. 

Consider a wave packet of size A propagating along the 
large-scale field B . We denote the corresponding per- 
turbations (i.e., typical variations across the eddy) of the 
velocity and magnetic fields by Sv\ and 8b\) in the Alfven 
wave, 8v\ ~ 5b\. Its interaction with the counter propa- 
gating packet of the same size occurs during time X/Va- 
As follows from (11121) . during one interaction the eddies 
arc deformed only slightly, ASv\ ~ {8v\/ X){X/Va)- Since 
different eddies are not correlated, the perturbations add 
up randomly, so the eddy is deformed considerably only 
after a large number of interactions, N ~ (5v/ ASv\) 2 . 
The time of energy transfer to a smaller eddy can thus be 
estimated as t ik (X) ~ NX/Va ~ \/5v\(Va/5v\). This 
time is larger than the Kolmogorov dynamic time, t(A) ~ 
\/6v\, by the Alfven factor Va/Sv\. Assuming that the 
energy flux over scales is constant, Sv^/tjk = const, we 
obtain the Iroshnikov-Kraichnan energy spectrum, 



E m (k) = (\5v(k)\ 2 )k 2 oc fc- 3/2 . 



(3) 



d t z + (V A ■ V)z + (w • V)z = —VP, 



(1) 



The essential assumption of the Iroshnikov-Kraichnan 
picture is that the eddy size is the same in the field- 



parallel and field-perpendicular directions. However, nu- 
merical and observational data accumulated for the last 
30 years indicate that in MHD turbulence the energy 
transfer occurs predominantly in the field-perpendicular 
direction, e.g., 0,0- This raises the question whether 
anisotropy is crucial for the energy cascade, and whether 
it changes the spectrum of turbulence. 

An elegant treatment of anisotropic MHD turbulence 
was proposed by Goldreich & Sridhar pj. They sug- 
gested that as the energy cascade proceeds to smaller 
scales, turbulent eddies progressively become elongated 
along the large-scale field. Their field-parallel and field- 
perpendicular scales are found from the so-called critical- 
balance condition. This condition follows from two dif- 
ferent estimates that are equivalent in the Goldreich- 
Sridhar picture. First, the field-parallel scale of an eddy 
is found from formal balance of the linear and nonlin- 
ear terms in the MHD equations 11121) . Va/1 ~ Sv\/X. 
Second, the field-parallel scale of an eddy can be ob- 
tained from the requirement that the magnetic field-line 
displacement in the eddy, £ ~ 5b\l/VA, be comparable 
with the field-perpendicular eddy size, A. The shape 
of the turbulent eddy in the Goldreich-Sridhar theory 
is schematically presented in Fig. As a result, two 
counter propagating eddies are deformed strongly dur- 
ing only one interaction, and the energy transfer time is 
given by the Alfven crossing time, tqs ~ I/Va ~ X/Sv\. 
The Goldreich-Sridhar theory thus predicts that due to 
local anisotropy, the energy-transfer time is reduced to 
the Kolmogorov estimate. The field-perpendicular en- 
ergy spectrum is obtained from the condition of constant 
energy flux, Sv 2 /tqs = const, which gives 

E G s(k±) = (\6v(k ± )\ 2 )k ± cx fc~ 5/3 , (4) 

where Sv(k±) = j 5v(xj_) exp(— zk^ ■ x±)d 2 x±. 

Recent high-resolution numerical simulations of MHD 
turbulence in a strong external magnetic field indeed con- 
firmed the elongation of turbulent fluctuations along the 
large-scale magnetic- field [fl 0, El • However, the field- 
perpendicular energy spectrum was consistently found 
to be close to E(k x ) oc fc~ 3/2 d H H E3- Obvi- 
ously, such a spectrum combined with the anisotropy of 
fluctuations contradicts both the Iroshnikov-Kraichnan 
and the Goldreich-Sridhar phenomenologies. This con- 
troversy motivated our interest in the problem. 

In this paper we argue that filament-like eddies are, 
in fact, non-realizable. We propose that the small-scale 
turbulent eddies spontaneously develop angular align- 
ment of their magnetic-field and velocity-field polariza- 
tions, which leads to their local anisotropy in the field- 
perpendicular plane. This effect is similar to the dynamic 
alignment known in the case of decaying MHD turbu- 
lence, where magnetic and velocity fluctuations approach 
the configuration v(x) = h(x ) or y(x) = — b(x), depend- 
ing of the initial conditions |14l Il5l llq . In the aligned 
state, the nonlinear interaction is zero, see Eqs. (11 121) . 
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FIG. 1: Sketch of a turbulent eddy in Goldreich-Sridhar pic- 
ture. The large-scale magnetic field is in the vertical direction. 
The field perpendicular dimensions of the eddy are the same, 
while its field-parallel scale is I oc A 2//3 . As the turbulent cas- 
cade proceeds toward the smallest, dissipative scales, A — > 0, 
the Goldriech- Sridhar eddy assumes the shape of a filament. 

We propose that in the case of driven turbulence the 
tendency to dynamic alignment is preserved, however, 
the precise alignment cannot be reached. The reason 
is an energy cascade toward small scales, which should 
be maintained by nonlinear interaction. We thus argue 
that at each scale A, the alignment of fluctuations should 
reach the maximal level consistent with a constant en- 
ergy flux through this scale. We demonstrate that this 
is achieved when the velocity and magnetic-field fluctu- 
ations 5v\ and ±5b> align their directions within the 
angle <j)\ cx A 1 / 4 . The dynamic alignment in driven tur- 
bulence thus becomes scale-dependent. Quite remark- 
ably, this leads to the field-perpendicular energy spec- 
trum E(k±) oc kj 3 ^ 2 , which explains the numerical ob- 
servations and resolves the above mentioned controversy. 

As another important result, in our theory small-scale 
eddies can be viewed as sheets or "ribbons", stretched 
along the magnetic-field lines. This explains the well 
known numerical fact that the dissipative structures 
in MHD turbulence are micro current sheets rather 
than filaments, e.g., 0, 0, In the next section we 
introduce our model of anisotropic MHD turbulence. 
Preliminary results on the dynamic alignment in driven 
MHD turbulence can be found in our earlier work |l7|. 

2. Structure and spectrum of MHD turbulence. — As 
one can check, the MHD equations HI2[) conserve the 
integrals J z 2 d 3 x and J w 2 d 3 x, if the fluctuations w(x) 
and z(x) have periodic boundary conditions or vanish at 
infinity. These integrals can be expressed through the 
integral of energy 

E=±J(b 2 + v 2 )d 3 x, (5) 

and the integral of cross-helicity, 

H c = /(v b)d 3 x. (6) 



In the unforced case, both integrals decay due to small 
viscosity and resistivity of the fluid. However, dissipa- 
tion of cross-helicity is not sign-definite, and, therefore, 
the integral of cross-helicity decays slower than the in- 
tegral of energy, e.g., 0. As a result of such "selective 
decay," turbulence approaches the perfectly aligned con- 
figuration h(x) = v(x) or h(x) = — v(x) depending on 
the initial conditions. This behavior is known as the dy- 
namic alignment or the Alfvenization effect [lj, lid llfi| . 
In the aligned state, either z(x) or w(x) is identically 
zero and nonlinear interaction vanishes. 

We propose that a similar effect is present in driven 
MHD turbulence, since the external force locally pro- 
duces large-scale fluctuations of cross-helicity, which are 
then inherited by smaller-scale eddies. Both E and H c 
cascade toward small scales, however, the cascade rate 
of cross-helicity may generally be smaller than that of 
energy, which forces fluid and magnetic fluctuations to 
align their polarizations at each given scale. However, 
the precise alignment cannot be reached, it would be in- 
consistent with the constant energy flux over scales. In- 
stead, the alignment of fluctuations should saturate at 
the maximal level that can be achieved in the presence 
of such a flux. 

Let us first describe the shape of the eddy, which would 
be dictated solely by a constant energy flux, without 
any constraints imposed by the cross-helicity conserva- 
tion (this derivation was first proposed in 17]). Assume 
that directions of shcar-Alfven velocity- and magnetic- 
field fluctuations 5w\ and ±<5b.\ are aligned within some 
(small) angle 9\ in the field-perpendicular plane. As one 
can directly check, this leads to depletion of the nonlin- 
ear interaction in Eqs. 11121) : (w • V)z ~ (z ■ V)w ~ 
9\5v\/\. Similarly to the Goldreich-Sridhar critical bal- 
ance, the eddy elongation in the field-parallel direction 
is found from balancing the linear and nonlinear terms 
in Eqs. I|1I2JI . / ~ Va\/(Sv\9\). The energy trans- 
fer time is then calculated as the Alfven crossing time, 
tn ~ I/Va ~ ^/{8v\8\). It is important that such tur- 
bulence is strong and essentially three-dimensional. 

To determine the shape of the eddy, we require that 
the energy flux be constant for all scales, 5v\/tn = const. 
This leads to the scaling of velocity fluctuations 8v\ cx 
(X/9\) 1 ^ 3 . The displacement of magnetic-field lines is 
given by £ ~ 8v\1/Va, and the correlation length of 
fluctuations in the field-displacement direction cannot be 
smaller than £. Remarkably, the obtained shape of the 
eddy satisfies A/£ ~ 9\, so it is indeed consistent with 
the assumed alignment of fluctuations within the an- 
gle 9\. Note that in contrast with the Goldreich-Sridhar 
picture, in our model the eddy is three-dimensionally 
anisotropic, I ^> £ ~> A, see Fig. 

It is natural to assume that turbulent fluctuations are 
scale invariant, which means that 9\ is a power-law func- 
tion of A. We may parametrize 9\ cx \ a ^ 3+a \ which 
leads to 6v x oc A 1 /^), £ oc A 3 /( 3+Q \ I cx A 2 /( 3+Q ). 




FIG. 2: Anisotropic turbulent eddy in our picture. The large- 
scale magnetic field is in the vertical direction. The field- 
perpendicular dimensions of the eddy are A and £ cx A 3,/4 , 
and the eddy size in the field-parallel direction is I oc A 1//2 . As 
the energy cascade proceeds toward the smallest, dissipative 
scales, A — » 0, the eddy assumes the shape of a current sheet. 



We thus obtain that the sole requirement of constant en- 
ergy flux does not define the eddy shape uniquely, but 
leads to a one-parameter family of solutions. The the- 
ory is self-consistent for an arbitrary parameter a > 0. 
(Note that Goldreich-Sridhar model is a particular so- 
lution corresponding to a = 0.) In order to address the 
crucial question about the value of a, we now have to use 
the second conserved quantity - cross-helicity. In other 
words, we want to find a that minimizes the total an- 
gular mismatch between the velocity and magnetic-field 
polarizations in the eddy. 

The mismatch angle in the field-perpendicular (hori- 
zontal) plane is 9\ cx \ a /^+ a \ However, the polariza- 
tion vectors are also mismatched in the vertical direc- 
tion. To obtain the vertical alignment angle, 9\, we note 
that in the regime of strong turbulence, eddies propagat- 
ing along a large-scale magnetic field interact efficiently 
during only one crossing time. Therefore, only the lo- 
cal direction of the magnetic field matters, and when we 
speak about eddy elongation in the field-parallel direc- 
tion, Z, we should mean the eddy dimension along the 
local magnetic field (this was established by Cho & Vish- 
niac |||). It is however important to note that the di- 
rection of the local magnetic field at the scale A can- 
not be defined precisely. Since the corresponding eddy 
contains magnetic field lines wandering within the angle 
9\ ~ £/l cx \ 1 /( 3+a \ the direction of the local magnetic 
field can only be defined with the same accuracy. This 
means that the directions of shear-Alfven velocity-field 
and magnetic-field fluctuations are aligned in the verti- 
cal direction within the angle 9\, as is sketched in Fig.|3 
Since both alignment angles, 9\ and 9\, are small, the 
total angular mismatch between 5v\ and ±<5bx can be 

calculated as <b\ - 



Following our strategy, we now require that the align- 
ment angle <fi\ be minimal. We however observe that 
the obtained shape of the eddy precludes us from achiev- 
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FIG. 3: Sketch of three-dimensional angular alignment of 
shear-Alfven velocity and magnetic-field fluctuations. The 
alignment angles consistent with an energy cascade are given 
by 9 cx A q/(3+q ° and oc A 1/(3+Ql) . The maximal alignment is 
achieved for a = 1 (see the text). 



ing the perfect alignment, <p\ = 0. Indeed, if for a given 
small scale A, we try to maximally align the polarizations 
in the field-perpendicular (horizontal) direction, i.e., to 
minimize 0\ oc X a /( 3+a ' : we need to set a — > oo. In 
this case, the fluctuations will be completely misaligned 
in the vertical direction, 6 oc A 1 ^ 3+Q ' ~ 1. Similarly, if 
we try to maximally align them in the vertical direction, 
a — ► 0, they become misaligned in the horizontal plane. 
This "uncertainty" is minimized when 6\ ~ 6\, in which 
case the maximal angular alignment is achieved and pre- 
served for all scales. This determines the scaling param- 
eter uniquely: a — 1. The resulting scaling of velocity 
fluctuations is 5v\ cx A 1 / 4 , and the field-perpendicular 
energy spectrum has the form 



E(k±) cx fc~ 3/2 . 



(7) 



The obtained structure and spectrum of turbulent 
fluctuations is the main result of this paper. 

3. Discussion and conclusion. — It may be reasonable 
to believe that an external magnetic field is not essen- 
tial for our derivation. Indeed, a local guiding field for 
small-scale fluctuations is naturally provided by large- 
scale eddies, e.g., 0,0- By this analogy, the spectrum 
of isotropic MHD turbulence should have scaling (JJJ 
as well. We however note that to observe this spec- 
trum in numerical simulations of isotropic turbulence one 
would need to reach extremely high resolution (to en- 
sure Sb\/Sb\ < 1), which is impossible with present-day 
computer power. 

We also note that our theory naturally explains the 
presence of ribbon-like dissipative structures (current 
sheets) in numerical simulations of MHD turbulence 0, 
8]. Indeed, the form of the eddy predicted in our model 
converges to such a structure as A — > 0. 

On the observational side, MHD turbulence is invoked 
to explain solar- wind measurements, e.g., and inter- 
stellar scintillation, e.g., [l9|. Although the inferred spec- 
tra of magnetic-field and electron-density fluctuations are 
broadly consistent with the —5/3 scaling, there do exist 
indications in favour of "—3/2" in some diffractive scin- 
tillation m. 



In conclusion, we propose that similarly to decaying 
MHD turbulence, driven MHD turbulence tends to align 
the polarizations of magnetic- and velocity-field fluctua- 
tions. However, the dynamic alignment cannot be per- 
fect: perfectly aligned fluctuations do not interact and 
cannot carry energy flux. We therefore require that the 
alignment be maximal under the constraint of constant 
energy flux. Such requirement defines the alignment an- 
gle uniquely, <f)\ oc A 1 / 4 , which means that the strength 
of nonlinear interaction in driven MHD turbulence is re- 
duced by the factor oc A 1 / 4 compared to a simple dimen- 
sional estimate (Sv\) 2 /X. The resulting fluctuations are 
three-dimensionally anisotropic (Fig.|2Jl, and their energy 



spectrum is E(k±) oc k 
merical results. 



-3/2 



in good agreement with nu- 



I am grateful to Peter Goldreich and Samuel Vainshtcin 
for important discussions. This work was supported by 
the NSF Center for Magnetic Self-Organization in Lab- 
oratory and Astrophysical Plasmas at the University of 
Chicago. 



[1 

[2 

[3. 
[4. 
[5 

[7; 

[8] 
[9] 

[10 

[11 

[12 

[13 
[14 

[is; 

[16 

[17 

[18 

[19 
[20 



Biskamp, D., 2003, Magnetohydrodynamic Turbulence. 
(Cambridge University Press, Cambridge). 
Iroshnikov, P. S., AZh, 40 (1963) 742; Sov. Astron, 7 
(1964) 566. 

Kraichnan, R. H., Phys. Fluids, 8 (1965) 1385. 
Goldreich, P. & Sridhar, S., ApJ 438 (1995) 763. 
Biskamp, D. & Muller, W.-C, Phys. Plasmas 7 (2000) 
4889. 

Cho, J. & Vishniac, E., ApJ. 539 (2000) 273. 
Milano, L. J., Matthaeus, W. H., Dmitruk, P., & Mont- 
gomery, D. C, Phys. Plasmas 8 (2001) 2673. 
Maron, J., & Goldreich, P., ApJ 554 (2001) 1175. 
Muller, W.-C, Biskamp, D., & Grappin, R., Phys. Rev. 
E 67 (2003) 066302. 

Muller, W.-C. & Grappin, R., Phys. Rev. Lett., 95 (2005) 
114502. 

Haugen, N. E. L., Brandenburg, A., & Dobler, W., Phys. 
Rev. E 70 (2004) 016308. 

Ng, C. S., Bhattacharjee, A., Germaschewski, K., & 
Galtier, S., Phys. Plasmas, 10 (2003) 1954. 
Galtie r, S., Pouquet, A., & Mangeney, A., 
physics/0504207 (2005). 

Dobrowolny, M., Mangeney, A, & Veltri, P., Phys. Rev. 
Lett. 45 (1980) 144. 

Grappin, R., Frisch, U, Leorat, J., & Pouquet, A., As- 
tron. Astrophys. 105 (1982) 6. 

Pouquet, A., Frisch, U, & Meneguzzi, M., Phys. Rev. 

A 33 (1986) 4266. 

Boldyrev, S., ApJ 626 (2005) L37. 

Goldstein, D. A. Roberts, D. A., & Matthaeus, W. A., 

Annu. Rev. Astron. Astrophys. 33 (1995) 283. 

Lithwick, Y. & Goldreich, P., ApJ 562 (2001) 279. 

Shishov, V. I., Smirnova, T. V., Sieber, W., Malofeev, V. 

M., Potapov, V. A., Stinebring, D., Kramer, M., Jessner, 

A., and Wielebinski, R., Astron. Astrophys. 404 (2003) 

557. 



